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Abstract. The identification of point sources poses a great challenge for the high 
energy community. We present a new approach to evaluate the likelihood of a set 
of sources being a Galactic population based on the simple assumption that galaxies 
similar to the Milky Way host comparable populations of gamma-ray emitters. We 
propose a luminosity constraint on Galactic source populations which complements 
existing approaches by constraining the abundance and spatial distribution of any 
objects of Galactic origin, rather than focusing on the properties of a specific candidate 
emitter. We use M31 as a proxy for the Milky Way, and demonstrate this technique by 
applying it to the unidentified EGRET sources. We find that it is highly improbable 
that the majority of the unidentified EGRET sources are members of a Galactic halo 
population (e.g., dark matter subhalos), but that current observations do not provide 
any constraints on all of these sources being Galactic objects if they reside entirely 
in the disk and bulge. Applying this method to upcoming observations by the Fermi 
Gamma-ray Space Telescope has the potential to exclude association of an even larger 
number of unidentified sources with any Galactic source class. 



1. Introduction 

The Energetic Gamma Ray Experiment Telescope (EGRET) measured gamma-ray 
emission at energies greater than 100 MeV across the entire sky. More than half of 
the sources in the third EGRET catalog []J were unidentified at the time of publication, 
and since then only a handful of those sources have been associated with known low- 
energy counterparts. Theoretical work has produced many candidate sources for these 
detections, including known Galactic source classes not previously confirmed as gamma- 
ray emitters such as microquasars Q , and newly proposed populations of high energy 
Galactic sources such as dark matter annihilation in subhalos [e.g., 
and in mini-spikes around intermediate mass black holes [lOj. 
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Figure 1. Spectral and variability indices and associated uncertainties for sources 
originally unidentified in the 3EG catalog. Spectral data is as published in the 3EG 
catalog; variability data is taken from Nolan et al. 



The most direct approach to the confident association of gamma-ray sources 
with suitable counterparts is multiwavelength follow-up of individual objects [e.g., 
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However, the sources in the third EGRET 
catalog typically have large error boxes, making association with known sources by 
positional coincidence difficult, and the large number of unidentified sources makes 
multiwavelength follow-up for every object impractical. 

Spectral and variability information can help by strengthening or excluding possible 
associations [e 
variability date 
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However, large uncertainties in EGRET's spectral and 
particularly for the faintest sources, hinder the use of this information 
to extract sub-populations within the unidentified sources [e.g., |23|, [24], |25|, |26j. Figured] 
shows the distribution of the initially unidentified EGRET sources in the parameters S 
(variability index from 27|) and a (spectral index from the 3EG catalog) with associated 
uncertainties. Although in principle association of sources with members of some source 
classes may be possible using correlations in these parameters, from Figure [D it is 
clear that any clustering in the a — 5 plane would be significantly smeared by the 
measurement uncertainties. For this reason, using such clustering to pick out distinct 
classes of gamma-ray emitters is very difficult with the EGRET data, although future 
experiments that can probe smaller variability timescales and better determine source 
spectra may be able to use this information to extract source populations. 

The situation is less favorable for associating sources with members of theoretically- 
motivated classes of gamma-ray emitters without previously established detections. 
The uncertainties inherent in the emission features of unconfirmed source classes 

§ Note that here we refer to general variability properties, such as variability indices, rather than 
correlated multiwavelength variability which, if present, can confirm the identification of a gamma-ray 
source with a lower-energy counterpart. 
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make association of individual sources tentative at best. Association of unidentified 
sources with members of new source classes by positional correlation with low-energy 
counterparts has been proposed [28|, [29(, but these techniques can produce spurious 



results 30]. Furthermore, approaches to associating sources with new populations 
which rely on matching unidentified gamma-ray sources with objects detected in 
other wavelengths cannot be used if the candidate counterparts are not necessarily 
detectable in other energy ranges (e.g., the case of annihilating dark matter clumps). 
While emission from dark matter annihilation may have unique signatures such as a 
"smoking gun" spectral line, for many scenarios it is unlikely that even next-generation 
experiments could detect such a signal. 

One way to make progress in making source associations that does not require any 
knowledge of candidate counterparts is the log N — log 5 diagram, the distribution of the 
number of sources with observed flux [e.g., |24|, |25j, |3jJ, |32|, |33j, |34j. However, variability, 
large uncertainties in measured fluxes, and the presence of multiple populations within 
the unidentified sources 32| pose problems for this method. Still, one advantage of 
using log N — log S distributions is that they enable comparison between Galactic and 
extragalactic populations without requiring a large number of detailed assumptions. 
While this is a population-level approach, rather than one which focuses on individual 
objects, it can help distinguish between Galactic and extragalactic populations, and 
thus provides information not only about the source populations themselves but also 
about possible contributions from unresolved members of the same source class to diffuse 
backgrounds [35 . 

The source distance is perhaps the most obvious difference between extragalactic 
and Galactic populations, and is an essential ingredient in connecting the source 
luminosity, which is often motivated by theoretical arguments but not directly 
observable, with the measured flux. Previous work has considered the expected distances 
and luminosities implied by associating unidentified sources with specific Galactic 
populations [25|, |36( . However, a challenge for this approach is finding an appropriate 
"standard candle" with which to compare the derived luminosities. In this work we 
revisit this type of approach, and introduce a novel standard candle with which we can 
compare the total luminosity of a proposed Galactic population: the total luminosity of 
M31. By assuming that the nearby galaxy M31 hosts gamma-ray emitting populations 
similar to those of the Milky Way, we determine whether the bulk of the unidentified 
sources in a given catalog can be of Galactic origin and place constraints on the spatial 
distribution of any proposed population without assuming a specific candidate emitter. 
Our approach uses only angular position and flux information to evaluate the plausibility 
of Galactic distributions of these sources. In anticipation of the large source catalogs 
expected from the recently launched Fermi Gamma-ray Space TelescopdjJ] and other 
experiments in the near future, we demonstrate the potential of this technique by 
applying it to the existing EGRET data set. 

|| http://fermi.gsfc.nasa.gov 
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Confident association of high-energy sources with counterparts is a difficult problem 
involving many degeneracies. As a result, individual constraints on the nature of 
detected sources are often weak. For this reason, as many independent constraints 
as possible are needed to robustly determine the nature of the unidentified sources. 
The approach presented here, which uses M31 as a standard candle for the gamma-ray 
luminosity of normal galaxies, involves its own uncertainties which we discuss in detail 
below. But, these uncertainties are independent of those of other methods, so in this 
way our proposed constraint adds a new piece to the puzzle. In $2] we discuss the set 
of unidentified sources from the third EGRET catalog, outline our approach for placing 
constraints on unidentified source populations, and present our results. We discuss the 
potential of this approach in the Fermi era in §|3] and conclude in §HJ 



2. Constraining a Galactic population with M31 

Like the Milky Way, M31 is a luminous, high surface brightness galaxy with a mass 
of ~ 10 12 M . Having formed and evolved in a similar environment, M31 is also akin 



to the Milky Way in its structure and dynamical properties [see, e.g., |37|]. Because of 
this close resemblance, it seems likely that any gamma-ray emitting source populations 
(including both known astrophysical source populations and proposed populations such 
as dark matter subhalos) found in the Milky Way are also present in M31. Under 
this assumption we test candidate Galactic populations by asking whether the total 
luminosity of a comparable population in M31 is consistent with current observational 
constraints on the total gamma-ray luminosity of M31. For simplicity, we assume that all 
point sources in M31 are at precisely the distance of M31, i.e., the collective luminosity 
of the point sources can be compared to the total luminosity of M31. 

EGRET did not detect M31, determining a 2a upper limit on the flux at the 
coordinates of M31 of ify,6) M31 < 4.9 x 1(T 8 cm" 2 s" 1 Blom et al. derive 
a more restrictive upper limit for M31 from archival EGRET data using a maximum 
likelihood technique and simultaneously fitting point source models of the M31 emission 
and nearby catalog sources. Their analysis finds a 2a upper limit for the emission from 
M31 of F M31 (>100 MeV) < 1.6 x KT 8 cm" 2 s _1 . Since this upper limit is roughly 3 
times smaller than the value given in Hartman et al. we conservatively adopt the 
limit from Hartman et al. []]]. We take the distance to M31 to be G?M3i = 785 kpc 39 



which gives the upper limit on the photon luminosity of M31 for energies greater than 
100 MeV, 

-^M3i,tot < 3.6 x 10 42 s™ 1 M31 total emission limit (1) 

or, in more convenient units, i^M3i,tot < 0.38 (kpc/cm) 2 s _1 . 

In determining the appropriate limit to place on the luminosity of candidate source 
populations, we note that the entire gamma-ray luminosity of M31 cannot be attributed 



to point sources. Pavlidou and Fields 4Qj calculate the expected flux of diffuse gamma- 
ray emission due to cosmic ray interactions to be Fdiff(>100 MeV) ~ 1.0 x 10~ 8 cm -2 
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s 1 , more than 20% of LM3i,tot- Consequently, the expected luminosity of M31 due to 
point sources is restricted to 

^M3i, P t < 2.9 x 10 42 s" 1 M31 point source limit (2) 

or -LM3i,pt ~ 0-30 (kpc/cm) 2 s" 1 . These upper limits constrain the gamma-ray emission 
of a point source within a solid angle set by EGRET's angular resolution. For E > 
100 MeV, the fraction of photons contained within a cone of half-angle 0psf = 3.3° is 
f 1 = 0.69, so the luminosity of a proposed population relevant for comparison with these 
upper limits is reduced by a factor of / 7 . We use these luminosity limits to constrain the 
Galactic distribution of the unidentified sources by requiring that the total luminosity 
of the unidentified sources L tot be consistent with these bounds, i.e., 

fjL tot < Lm3i (3) 

taking Lm3i to be one of the upper limits defined in Equations CD and [2j We note that 
^M3i,pt is dependent on the adopted model for the truly diffuse emission, but we include 
it for completeness. 

Emission from unresolved point sources is another component of the Galactic 
luminosity which is measured as diffuse emission. The contribution from unresolved 
sources may be relevant because our assumption that M31 and the Milky Way host 
similar populations of gamma-ray emitting sources requires that unresolved sources 
enhancing the Galactic diffuse emission also be present in M31. The prediction we use for 
the diffuse flux of M31 only represents the expected genuinely diffuse emission, and does 
not account for the luminosity enhancement by unresolved point sources. In particular, 
it is likely that unresolved members of the same class or classes as the unidentified 
sources are present and contribute to the diffuse luminosity of the Milky Way, if these 
populations are indeed Galactic. However, in order to appropriately account for this 
emission, it would be necessary to make numerous assumptions about the properties of 
the proposed Galactic population, including the cumulative flux distribution extending 
below the flux sensitivity of EGRET and the location of the unresolved sources. 
Without assuming a specific candidate emitter, we have no theoretical basis on which to 
construct either of these inputs. Although one might attempt to determine the angular 
distribution of the unresolved sources based on observations of the Galactic diffuse flux, 
the distances (and hence the collective luminosity) of the unresolved sources would be 
highly model-dependent and sensitive to assumptions about the spatial distribution of 
the genuinely diffuse emission. Consequently, due to large uncertainties in the magnitude 
of this emission, we again make the very conservative choice to ignore this guaranteed 
component when placing constraints, allowing all of the expected point source luminosity 
to be attributed to the resolved unidentified sources. 

Additionally, there is already a guaranteed contribution to the luminosity of point 
sources from confirmed detections of members of known Galactic source classes, such as 
pulsars and supernova remnants, which further reduces the luminosity available to the 
unidentified sources. However, the total gamma-ray luminosity of these objects is very 
small compared to the diffuse luminosity of either the Milky Way or M31. Including 
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this component would not alter our results in any substantial way, so for simplicity we 
ignore it. 

In this work we ask whether the unidentified sources (or a subset of them) could 
likely be members of a Galactic population. Since the flux, rather than the luminosity, 
of the unidentified sources is known, it is necessary to assume a distance to each source 
to determine the luminosity of a proposed population consisting of these sources. Most 
candidate Galactic populations are likely to be associated with either the disk and bulge 
(e.g., pulsars, supernova remnants) or the dark matter halo (e.g., dark matter clumps, 
intermediate mass black holes), so in these cases we have a natural expectation for the 
spatial distribution of the sources. In particular, for the case of the disk and bulge 
populations, we expect that the typical distances of the sources will reflect the spatial 
extent of the disk, roughly 30 kpc in radius (as in the emission model of 4lj )• For 
populations associated with the dark matter halo, we anticipate a distribution of the 
source distances out to at least the virial radius of the Galaxy r vir , which is 258 kpc 
in the favored model of Klypin et al. [37], and note that recent numerical simulations 
find that the dark matter subhalo distribution extends significantly beyond this distance 



42j, |43j, M . 

The angular distribution of sources can be an important discriminator between 
various Galactic distributions of sources and extragalactic sources. Due to the relatively 
small offset of our position from the Galactic Center compared to the typical distances 
of sources associated with the dark matter halo, the angular distribution of such sources 
would appear largely isotropic. In contrast, sources associated with the disk and bulge 
would be concentrated at low latitude. However, due to EGRET's unequal exposure 
map, using isotropy information to constrain the EGRET unidentified sources requires 
careful study and is beyond the scope of this work. 



2.1. The unidentified sources 

The third EGRET catalog contains 271 sources (E>100 MeV), 101 of which were 
initially identified or associated with likely lower-energy counterparts by the EGRET 
team. The majority of these sources are blazars, along with a smaller number of 
pulsars, the Large Magellanic Cloud, a radio galaxy, and a solar flare. An additional 
48 of the originally unidentified sources subsequently have been suggested to be 
associated with plausible counterparts, including binaries, supernova remnants, gas 
clouds, microquasars, and black holes. We evaluate the total luminosity of possible 
Galactic distributions of the unidentified sources using two subsets of sources from the 
EGRET catalog as our candidate populations, and emphasize that the hypothesis we 
are attempting to reject with the luminosity constraint is that all of the unidentified 
sources are of Galactic origin. 

We first consider the set of all sources published without identifications in the 
catalog. To avoid biasing our sample toward higher fluxes, for each source we use the 
P1234 flux (summed flux over cycles 1, 2, 3, and 4). As a result, we exclude two 
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additional sources which were not published with P1234 fluxes, leaving a total of 168 
sources in this set, which we term the "original sample". 

We also construct a more restricted source sample by excluding those sources for 
which plausible associations with low energy counterparts have been suggested since 
the catalog was published. We make use of an updated listing of EGRET unidentified 
sources, and exclude from our sample all sources for which an association has been 
suggested, regardless of the significance of the association. Associations included in 
this compilation are simply the results reported in recent publications; the validity of 
these suggested counterparts has not been evaluated by a single standard. We also 
remove from this sample the few sources noted by the EGRET team as possible or 
likely artifacts. After applying these restrictions, our "restricted sample" consists of 
119 unidentified sources. Omitting these sources from our restricted sample is again a 
conservative choice because increasing the number of sources in our candidate population 
would lead to an increase in the total population luminosity, and hence any distribution 
would be more likely to violate the M31 luminosity constraint. We comment that 



Sowards-Emmerd et al. [45|, |46| report plausible associations of many high-latitude 
unidentified EGRET sources with blazars, and these associations are not reflected in 
our restricted sample. We also note that a revised EGRET source catalog has been 
published by Casandjian and Grenier 47J that does not confirm the detections of several 
3EG sources. Although Fermi has now released a bright gamma-ray source list [48], we 
do not consider those sources in this study since they are high-confidence, high-flux 
sources, and most of the Fermi unidentified sources will likely be faint, resulting in 
larger positional error circles. 



2.2. Simple luminosity tests 

The EGRET data provides a measurement of the flux and the angular position of 
each source, so the distance to each unidentified source is needed to calculate the total 
luminosity of a candidate population: 

i i 

where the summation is over each unidentified source i and the flux, luminosity, and 
distance of each source are denoted by F iy L iy and di respectively. We begin by using 
some simple models to assign source distances. These models are not realistic, but 
instead indicate the typical allowed distances of these objects from our position in the 
most extreme scenarios. 

We first consider the case that all unidentified sources are at a fixed distance d 
from us, i.e., di = do, and use Equation @] to determine the total population luminosity 
L tot in terms of d . We then ask what is the maximum distance do at which we can 
place the sources without violating the M31 luminosity limit IvM3i,tot? Because sources 
observed to have a given flux will be more luminous if located at a greater distance from 
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the observer, the luminosity constraint f 1 L tot < LM3i,tot limits the maximum distance: 

i^i ^M31,tot / r x 

Using the original source sample we find that the maximum distance do = 30 kpc 
which is plausible for a disk and bulge source population, however, a population of halo 
sources would extend at least as far as r vir = 258 kpc, and consequently would have 
typical distances from our position much greater than 30 kpc. This simple test implies 
that the unidentified sources are not primarily associated with the dark matter halo of 
the Galaxy. 

We next consider a population of sources which all have the same intrinsic 
luminosity Lq. We consider the original sample as above to set the individual source 
luminosity L = L tot /N, with iV = 168 sources, and set the total relevant population 
luminosity f-yL tot equal to the total luminosity upper limit -^M3i,tot- Using the measured 
flux of each source, we calculate the corresponding distance to each source in this 
scenario, 

d i = .A5I= i /i^L. (6) 

Under this assumption the sources are at a mean distance of 37 kpc and we find, as 
for the fixed distance case presented above, that the typical distances of the sources are 
uncomfortably small for a plausible halo population, but quite reasonable for a disk and 
bulge population. 

Naturally, we now ask what is the maximum average distance d = Yli di/N allowed 
while still respecting the luminosity constraint? We extremize d subject to the condition 
that our set of fluxes Fi and assigned distances d-i produce a given total population 
luminosity L tot via Equation H] (which we will take to be Lu3i,tot / f-y) , arid find that the 
maximum average distance is given by 




For the original source sample and the total luminosity upper limit L]yi3i,tot, ^max = 40 
kpc. Once again, this distance is too small to be characteristic of a halo population, but 
consistent with a disk and bulge distribution. 

2.3. Assigning source distances 

We now take a more sophisticated approach to construct our population's spatial 
distribution. Our method for assigning distances is motivated by the goal of testing 
proposed source populations, so here we determine the distance to each source by 
considering the expected Galactic distribution of various types of populations. We 
emphasize, however, that alternative methods of assigning source distances can also be 
used to determine the luminosity of the candidate population, so the luminosity test we 
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propose can be used in conjunction with existing approaches to provide an additional, 
complementary constraint. In general, candidate Galactic populations can be classified 
as living either in the halo or in the disk and bulge. Because these populations can be 
associated with a measured mass component of the Galaxy, we use mass density as a 
proxy for source density. 

For source populations expected to reside in the disk and bulge, we approximate 



the mass distribution with a Miyamoto-Nagai disk 49j and a Hernquist bulge 50] . We 
take the disk mass to be Md = 4 x 10 10 M Q 371 ]. and the scale parameters a and b to 



be 6.5 kpc and 0.26 kpc respectively 51J. For the bulge mass we use Mb = 8.0 x 10 



M Q (mi + m 2 in [37]]), with the scale parameter a = 0.7 kpc [5JJ. The disk and bulge 



profiles are truncated at a radius of 30 kpc from the Galactic Center. 

For the case of candidate populations correlated with the dark matter distribution 
of the Galaxy, we model the dark matter halo of the Milky Way using a NFW profile 



52j | with mass Mh a io — 10 12 M and concentration c = 12 37j, truncating the profile 



at r v ; r = 258 kpc from the Galactic Center. We note that since r V i T represents a sizable 
fraction of the distance between the Milky Way and M31, the difference between the flux 
measured at the Milky Way from a source r v i r closer or farther from us than the distance 
of M31 is not negligible. In other words, the extent of M31's source population along 
the line-of-sight could be large enough that our assumption that all sources associated 
with M31 are at precisely the distance of the center of M31 is a crude approximation. 
The largest variations possible due to this effect are, however, small compared to the 
uncertainties inherent in our approach: in the most extreme cases an individual source's 
flux may be enhanced by up to a factor of ~ 2.2 or decreased by a factor of ~ 1.8. For 
an isotropically distributed population, the overall change in the M31 luminosity is not 
significant. 

For M31, the angle Opsf of 3.3° corresponds to a projected radius of ~ 45 kpc, 
which encloses the expected emission region of the disk and bulge, but does not fully 
enclose the dark matter halo. For the halo parameters used here, a cylinder with an 
angular radius of 3.3° centered on M31 encloses a fraction /halo ~ 0.38 of the halo 
mass, so for a halo distribution we assume the total source luminosity within the cone 
reflecting EGRET's angular resolution will be further reduced by a factor /halo beyond 
the factor / 7 reflecting the photon containment fraction. 

We use the P1234 fluxes and the angular positions of the sources as measured 
by EGRET. Each source in the catalog is reported with an estimate of positional 
uncertainty, 695, typically ~1°. This value is for most objects the angular radius of 
a circle containing the same solid angle as the 95% confidence level contour. However, 
in general the contour is not a circle but rather a complex shape, and so for simplicity 
we will consider B95 to be the angular uncertainty of both the Galactic longitude I 
and latitude b measurements. We note that Mattox et al. [53( have provided elliptical 
contour fits to the EGRET sources, but because we do not expect the mass density to 
vary significantly over such small scales, generalizing the error contours to squares is a 
suitable approximation for our purposes. 



Unidentified high energy sources 



10 




Figure 2. We calculate the cumulative distribution function for the distance to each 
object by integrating the mass enclosed within the error boxes from the observer along 
the line of sight. 

For each angular source position (£, b) we construct a cumulative distribution 
function for a given Galactic mass distribution to describe the likelihood of the source 
being within a particular distance along that line of sight. The probability of a source 
with a given angular position being located within a distance d from us is then given 
by V(d) = M eac i(d) / M enc i(d max ) , where d max is the distance along that line of sight at 
which we truncate our mass distribution to produce a finite volume in which to place 
the sources. The mass enclosed along the line of sight within a solid angle defined by 
the positional error boxes and extending from the observer out to a distance d is given 
by 

M end (d)= [ [ + [ + p(z,£,b)z 2 cos(b)dbd£dz, (8) 

Jo J t- Jb^ 

where £± = £ ± 695, b± = b ± ©95, and p is the mass density. The variable z describes 
integration along the line of sight. This is illustrated schematically in Figure [2j 

We use a Monte Carlo algorithm to generate realizations of the population 
distributions. For each realization a distance is assigned to each source by sampling 
the appropriate cumulative distribution function associated with the angular position 
of that source. The total luminosity for the realization is then calculated. 

2.4- Results 

Figure [3] shows the results of the luminosity test for disk and bulge distributions. For 
each of our samples, we first construct our realizations using only the sources which 
roughly appear to be located in the Galactic plane (|6| < 5°), and then using all sources. 
Fitting these distributions to a Gaussian we find the mean of the distribution \i = 0.040 
for the Galactic plane sources of our original sample, and \i = 0.053 for all sources in 
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Figure 3. Probability density functions for total luminosity of a population of sources 
associated with the disk and bulge. The distribution of total luminosities for 5000 
realizations, using sources with \b\ < 5° (light gray) and using all sources (black), is 
shown for the original sample set (left panel) and the restricted sample set (right panel). 
The M31 total luminosity upper limit £M3i,tot (dashed) and derived upper limit for 
the luminosity of point sources £M3i, P t (dotted) are shown for reference. 




Figure 4. Probability density function for total luminosity of a halo population of 
sources. The resulting distribution for 5000 realizations using only sources with \b\ > 5° 
(light gray) is shown along with the corresponding distribution when all sources are 
used (black), drawing from the original sample (left panel) and the restricted sample 
(right panel). The M31 limit on the total gamma-ray luminosity LM3i,tot (dashed) and 
derived limit for the luminosity of point sources £M3i, P t (dotted) are also shown. 



that sample, in units of (kpc/cm) 2 s _1 . For the restricted sample, the mean values for 
the Galactic plane sources and for all sources are \i = 0.014 and \i = 0.024 respectively. 
In all four cases L M 3i,tot an d L M31tPt lie far above the total luminosity distributions, thus 
our luminosity test is consistent with the possibility that all of the unidentified sources 
are part of a Galactic population distributed in the disk and bulge. We note that this 
finding is consistent with that of Grenier 54| , in which a complementary approach based 
on isotropy arguments was taken to determine whether the unidentified sources could 
be associated with various Galactic structures. 

Figure H] shows distributions of the total source luminosity for a halo population. 
Following a similar approach to that for the disk and bulge distributions, for this scenario 
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we evaluate the total population luminosity for both of our source samples by first 
restricting each sample to sources with \b\ > 5° (omitting sources which appear to 
be in the Galactic plane), and then considering all sources in each sample. Assuming a 
Gaussian distribution, the mean value for the original sample using sources with |fe| > 5° 
is ii — 0.64, in units of (kpc/cm) 2 s _1 . Without the latitude restriction, the distribution 
shifts up to fi = 1.7. Our more conservative, restricted sample also exceeds the M31 
limits in almost all realizations, with fi = 0.49 for sources with \b\ > 5°, and \x = 0.84 for 
all sources in this sample. Based on these results, the possibility of any halo distribution 
of these sets of unidentified sources producing a luminosity below even the maximum 
observational limit is strongly disfavored. 



3. Discussion 



With the launch of Fermi last year, this is an ideal time to consider source identification 
strategies. The number of detected gamma-ray sources will be increased substantially 
by Fermi, and with it the number of sources to be identified with known low-energy 
counterparts. With Fermi's expanded source catalog, following each detection with 
multiwavelength observations will certainly be an impractical method for making 
identifications. Although other gamma-ray telescopes currently in operation such as 
HESS, MAGIC, and VERITAS will be able to use their superior angular resolution 
and variability measurements to help identify sources [e.g.,|20j], they face similar issues 
with regard to the burden of multiwavelength follow-up observations of large numbers of 
objects not easily identified by other methods. In this work we proposed an additional 
constraint to aid in identifying populations of sources which is complementary to 
constraints from existing approaches. 



Of particular interest for this study is the likely detection of M31 by Fermi 40j, [55 



A measurement of the M31 gamma-ray flux will significantly narrow the range of 
plausible point source luminosities and provide stronger constraints on the properties 
of new Galactic gamma-ray emitting populations. Furthermore, Fermi may be able to 
provide detailed spatial and spectral information about the gamma-ray emission of the 
Milky Way and M31, which would test the assumption of similarity underlying this 
work. A determination of the M31 flux would also test the validity of the predicted 
diffuse gamma-ray emission. 

As an independent test of the plausibility of a candidate Galactic population, the 
method proposed here can be used in conjunction with isotropy studies to generate 
stronger constraints than can be achieved by isotropy arguments alone. Similarly, the 
luminosity constraint outlined here stands to benefit greatly from isotropy information, 
since satisfying the luminosity constraint alone does not guarantee that a proposed 
distribution is plausible. Any realistic Galactic population is expected to be distributed 
according to symmetries of the Galaxy (e.g., it is unlikely that a population would 
cluster near our position). Our method of assigning distances effectively projects all 
of the sources into the disk and bulge or into the halo, so additional constraints, such 
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as those from isotropy studies, are necessary to determine whether realizations of a 
population represent plausible spatial distributions. We have shown here that most of 
the EGRET unidentified sources cannot be distributed in the halo if the M31 constraint 
is to be satisfied, so this is a particularly important issue for sources with high Galactic 
latitude which, naively, are least likely to be associated with the disk and bulge. In 
general, as evident from Figures [3] and HI this method is intrinsically weaker for sources 
associated with the disk and bulge, which are typically much nearer to our position than 
halo sources. 

It is important to keep in mind that although we expect M31 and the Milky Way 
to be similar in gamma-rays, clearly we do not expect them to be identical. Depending 
on the characteristics of the source population and the properties of M31, we expect 
that the total population luminosities we calculated by placing the unidentified sources 
in the Milky Way could vary by a factor of a few for a corresponding population in M31. 

There is currently only an upper limit to the M31 flux, so a flux measurement 
by Fermi would necessarily affirm or tighten the luminosity constraint. Once a 
determination of the M31 flux is made, it is possible that no new Galactic populations 
could be accommodated, unless M31 and the Milky Way are fundamentally different 
in their gamma-ray emission properties. Furthermore, a flux determination would 
also suggest a lower limit to the point source flux. In this study, we did not take 
into account a number of probable contributors to the Galactic luminosity: diffuse 
emission from unresolved point sources, identified Milky Way objects, and unidentified 
sources with recently suggested Galactic counterparts. We expect that M31 will have 
similar populations of these Galactic sources, so in principle this will decrease the 
allowed luminosity of new proposed populations, and consequently further constrain 
their abundance and distribution. 

4. Conclusions 

We outlined a new approach to constraining the origin of unidentified gamma-ray 
sources based on a luminosity constraint. This approach provides a constraint on 
candidate populations which is complementary to, and can be used in conjunction with, 
constraints from existing approaches. We demonstrated this method by applying it to 
the unidentified sources in the third EGRET catalog. Using the observational upper 
limit for the gamma-ray luminosity of M31 along with the assumption that M31 and the 
Milky Way host similar gamma-ray emitting populations, we constrained the allowed 
spatial distribution of these sources. 

We found that it is highly unlikely that a substantial fraction of the EGRET 
unidentified sources are members of a Galactic halo population. For the case of source 
populations expected to be associated with the disk and bulge, we find that within 
reasonable uncertainties all of the EGRET unidentified sources can be of Galactic origin 
without exceeding our assumed luminosity upper limit, although isotropy studies could 
likely further constrain this scenario. 
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The new information Fermi brings, with respect to both known gamma-ray sources 
and possible detections of proposed emitters, will help to determine the nature of the 
unidentified sources. Statistical techniques such as those employed in this study will 
be particularly useful in making meaningful statements about gamma-ray emitting 
populations when large numbers of individual identifications are not feasible. 
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